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ABSTRACT 


A transient ballistics and combustion model is derived to represent the 
closed vessel experiment that is widely used to characterize propellants. The 
model Incorporates the nitramine combustion mechanisms which are contained within 
the steady-state model developed as a part of last year's (FY 1976) research. A 
computer program is developed to solve the time-dependent equations, and is 
applied to explain aspects of closed vessel behavior. It is found that the rate 
of pressurization in the closed vessel is insufficient at pressures of interest 
to augment the burning rate by time-dependent processes. In the case of nitramine 
propellants, however, the cratering of the burning surface associated with combustion 
above break-point pressures augments the effective burning rate as deduced from 
the closed vessel experiment. Low pressure combustion is significantly affected 
by the ignition process and, in the case of nitramine propellants, by the developing 
and changing surface structure. Thus, burning rates deduced from the closed vessel 
experiment may or may not agree with those measured in the equilibrium strand 
burner. Application of the model to closed vessel test cases shows good agreement 
between theory and experiment. 

Series of T-burner experiments are performed to compare the combustion instability 
characteristics of nitramine (HMX) - containing propellants and ammonium perchlorate 
(AP)propellants. It is found that the inclusion of HMX consistently renders the 
propellant more stable. Although ash produced by more fuel-rich propellants could 
have provided mechanical suppresssion, results from clean-burning propellants permit 
the conclusion that WiX reduces the acoustic driving. 

Additional strand burning rate data are presented which provide further confirma- 
tion of the results of the steady-state modeling effect of FY 1976. 
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SECTION 1 
OBJECTIVE 

The objective of this research program is to study two aspects of the transient 
combustion properties of nitramine smokeless propellants. One is to detemtine the 
role of the nitramine ingredient in the combustion driving of acoustic instability 
in solid rocket motors. The second is to determine the effect of a transient 
pressure environment on the combustion of nitramine propellants in armament or 
closed vessel apparatus. The steady-state combustion characteristics of nitramine 
propellants have been described as a result of work performed under AFOSR Support 
Agreements AFOSR-ISSA-75-0005 and AFOSR-ISSA-76-0006 (Ref. 1). The present research 
begins to investigate unsteady-state or transient combustion characteristics. 
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SECTION 2 
INTRODUCTION 

• 

Nitramlnes are of Interest for gun and rocket propellant applications because 
they are an energetic source of smokeless combustion products (Refs. 2-4). Recent 
studies have been devoted to understanding and Improving the steady-sitate combustion 
characteristics of nitramlne propellants (Refs. 1, 5-10). Rocket motors, are. 
generally, principally concerned with steady-state performance. Although It Is 
recognized that gun propellants perform In a highly transient pressure environment, 
there Is basis for the Inference that Improvements In steady-state combustion will 
be reflected as Improvements In the transient ballistics performance (Ref. 11). 

Anomalous combustion behavior can be brought about by unsteady processes occur- 
ring In combustion chambers. In rocket motors, the most prevalent behavior Is com- 
bustion Instability. Combustion Instability Is of particular concern In smokeless • 
propellants because the absence of aluminum powder deprives the system of significant 
stabilization (Ref. 12). Virtually no work has been done to describe the Influence of 
nitramlnes on combustion Instability properties. In projectiles, the extreme trans- 
ient pressure environment can alter the combustion process from the steady-state as 
a result of time lags In the combustion response. Although this problem has been 
studied extensively with respect to rocket motor Ignition and controlled termina- 
tion (e.g.. Refs. 13-17), It has not been applied to gun propellants or to the more 
extreme projectile environment. Burning rates derived from closed-vessel (transient 
pressure) testing differ from those measured at constant pressure, which raises a 
question about the predictability of ballistics performance (Ref. 18). Further, the 
transient pressure environment may (Ref. 11) or may not (Ref. 19) aggravate any 
steady-state combustion anomaly to Impair the quality of the Impetus delivered. 

It Is considered that understanding the transient combustion characteristics of 
nitramlne smokeless propellants will provide additional direction for future develop- 
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mental efforts. The mechanistic understanding derived from the steady-state research 
provides an essential background, and an understanding of the transient phenomena 
will afford improved tailoring and predictive capabilities. 
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SECTION 3 

TRANSIENT' COMBUSTION MODEL DEVELOPMENT 

3.1 REVIEW OF CLOSED VESSEL BURNING RATE DATA 

A review was conducted of available data wherein direct comparisons were made 
between burning rates as measured in a closed vessel (transient pressure) and as 
measured in a strand burner (equilibrium pressure) (Refs. 20-24). Some of the closed 
vessel data were accompanied by the associated test pressure-time, or pressure-dp/dt 
(P-dot) histories. As will be discussed below, it is important to know the P-dot 
values for transient analysis. Some of the data were for nitramine propellants; 
some were for homogeneous (unfilled single-base or double-base) active binder pro- 
pellants. In general, there appears to be very good agreement between the closed 
vessel and strand data for the homogeneous propellants. Such was the conclusion of 
Grollman (Ref. 23) in his study. Figures contained in Ref. (24) show a virtual 
identity between closed vessel and strand data. Lenschitz' data for homogeneous 
propellants (Ref. 20), some of which are reproduced here as Fig. 1, also show good 
agreement between the two methods. On the other hand, there is a significant and 
consistent difference between closed vessel and strand burning rates in the case of 
propellants containing HMX (Refs. 20-22). Examples from Ref. (20) are presented 
here as Fig. 2. At low pressure, the closed vessel rate tends to lag the strand 
rate; at high pressure, the closed vessel rate tends to lead the strand rate. As 
a result, the closed vessel data exhibit a more uniform and higher apparent pressure 
exponent. 


3.2 APPLICATION OF SIMPLIFIED "P-DOT” THEORY 

Where there are differences between closed vessel and strand data, such 
differences have been attributed to the deductive nature of the closed vessel 
experiment. The closed vessel experiment measures a variable pressure versus 
time whereas the strand experiment measures directly the distance burned versus 
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Comparison of Strand and Closed Vessel Data for 
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Comparison of Strand and Closed Vessel Data for 
Nitramine Propellants 
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time at a constant pressure. From the pressure versus time, ballistics analysis 
yields mass generation rate versus time If certain assumptions are made. If the 
geometry of the sample Is known, the burning rate versus time may then be calculated. 
Thus, burning rate may be cross-plotted versus pressure. However, It has been 
suggested that the assumptions may be Incorrect and that the geometry of the sample 
may not be known. Unknown factors frequently mentioned Include heat losses, non- 
uniform Ignition and burning, and charge break-up (Ref. 18). Low pressure uncertainties 
have limited the reporting of data to pressures above 5000 psi. 

A mechanistic basis for a difference between burning rates measured under 
steady-state and transient pressure conditions was developed In the course of 
studies of Ignition and depressurization transients In solid rocket motors (Ref. 13). 
This mechanism was the thermal wave relaxation time In response to a pressure transient. 
A simplified form of expression was as follows (Ref. 14): 



1 + 



P 

P 


r^ = Instantaneous burn rate at pressure, P 

r = steady-state burn rate at P 

n = steady-state pressure exponent at P 

K » propellant thermal diffusivlty 
• 

P = Instantaneous rate of change of pressure at P 


( 1 ) 


More detailed analysis resulted In the development of computer programs such 
as described In Ref. (16). Some of these programs have been reviewed by Kooker 
(Ref. 25) In the context of gun ballistics. 
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An order-of -magnitude analysis of Eq. (1) reveals that P-dot effects should 
be small In the usual closed vessel experiment. Typical values of the parameters 
are ad follows: 


r * 2 In. /sec 8 10.000 psi (Fig. 1) 
n ' 0.8 (Fig. 1) 

K *0.0002 In. ^/sec. (typical) 

P = 1000 P (Ref. 20) 

The Impact of P-dot upon burning rate falls below 10% at pressures above 
10 Kpsi, and decreases with Increasing pressures and burning rates. This result 
Is consistent with the data In Fig. 1. Although the values of P-dot are very 
high when compared to values encountered In rocket motors, the values of pres- 
sure and r-squared also become very high so as to diminish the effect. Since 
most of the reported closed vessel data are at such high pressures, It would 
appear that agreement between the closed vessel and strand data should be 
expected assuming that P-dot would be the only reason for any difference. This 
appears to be borne out by the homogeneous propellants. 

Results of P-dot calculations applied to the nitramlne propellant data are 
Included as the dashed lines In Fig. 2. It was possible to do so because the 
P-dot histories for these tests were available (Ref. 20). Note that the calcu- 
lated effect has the tendency to reduce the pressure exponent In the high slope 
region below 5 Kpsi. The calculated effect becomes small at high pressure, as 
discussed previously. However, the data behave quite differently from the cal- 
culated effect. In view of the fact that experimental uncertainties and P-dot 
effects appear to be adequately resolved In the case of homogeneous propellants. 
It follows that something different must be happening In the case of nitramlne 
propellants. 
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3.3 THE ANALYTICAL MODEL 

3.3.1 Conceptual Basis 

The analytical model for the transient combustion of nitramlne propellants 
Is based upon the nitramlne surface structure mechanisms contained within the 
steady-state model (Refs. 1. 7), coupled to a thermal wave relaxation model 
(Ref. 16) of a homogeneous solid. 

At low pressure, and following Ignition, a planar melt layer must form 
and develop. This development will be an endothermic contribution to the 
energy balance, will tend to oppose the ordinary collapsing of the thermal wave 
during a pressure rise, and will tend to delay a break point. These effects 
should produce a burning rate lag at low pressure. At high pressure, and where 
coarse particles are present, the development of the surface craters In accordance 
with break points will increase the burning surface area and thereby augment the 
mass burning rate. As discussed earlier, the closed vessel deduces the mass 
burning rate and does not measure directly the linear burning rate. Therefore, 
if the burning surface Is underestimated, the linear rate as computed from the 
closed vessel data will appear to lead. 

Examples of cratered surfaces appear In Fig. 9 of Ref. (1). If the cratering 
mechanism increases both the linear burning rate and the surface area of the 
small pellets used In transient devices, then the mass burning rate would be 
proportional to the square of this process. Thus, the exponent break point 
would have a more serious effect In pressurizing chambers. 

3.3.2 Basic Equations for a Homogeneous Propellant 

Instantaneous burning rate Is related to a reference or steady-state value 
by means of an Arrhenius equation for surface decomposition: 
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- r exp C(-E 5 /Rq)( 1 /T^ - l/Tg) ] (2) 

Ej » activation energy of surface decomposition 
Rq « universal gas constant 

= Instantaneous surface temperature 
Tg » surface temperature corresponding to r 
r » steady-state burning rate at P, known from data or 
the steady-state model 


The Instantaneous surface temperature Is calculated from the time-dependent 
Fourier equation: 


II 

3 t 





exp (-E^/RJ) 


(3) 


T = temperature distribution 
t = time 

X » distance into the solid 
Qss = condensed phase heat of decomposition 
A * prefactor 

c = propellant heat capacity 


For numerical solution, Eq. (3) Is written In a suitable finite difference 
form such that the nominal mesh size Is selectable depending upon the thermal 
wave thickness. The selection affects accuracy, stability and computer time. 
The Initial and boundary conditions for Eq. (3), using Ref. (1) as the basis 
for an energy balance, are: 
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T (x.O) = 

fl (-.t) . 0 


Qs - k S (O.t) - fe,t * 8''1 ['<n - ^o> * Oss] ' 




w = web thickness 

= initial bulk temperature 

f . = external heat source, for example as would be 
input to produce ignition 

k = propellant thermal conductivity 
« 

= heat flux at the propellant surface 
p = propellant density 

I = dimensionless flame height corresponding to r at P 
The pressure variations are computed by a mass and energy balance for the 




"c dt 


Tg » gas temperature 
m^ = mass of gas stored in chamber 
R » gas constant 
V = chamber gas volume 
n = covolume 


dT = "i 3t “ Vi 

Ajj ® burning surface area 


For a closed vessel, it is assumed that there is no mass exhaust. 
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dT_ dm^ 

"c w ■ - ^g> -ar • 

Y " ratio of specific heats 

Tp * adiabatic flame temperature (a function of P) 

Cg » gas heat capacity (a function of P) 

0^ = heat loss rate to wall 

V\, * '■cfTg - T«) * (-), (t/ - t/) (8) 

= wall area 

h^ » suitable heat transfer coefficient 
» wall temperature 

(eo)g = product of gas emissivlty and Stefan - Boltzmann constant 

The unsteady heat transfer In the wall Is solved In analogous fashion to 
Eqs. (3) and (4), using Eq. (8) to describe the flux at the boundary; but r^=0 and 
Qss*0 ^or the wall. 

The foregoing equations, plus Eq. (3) written for the wall, constitute a system 

• • 

for the unknowns P, V, m^, Tg, Q^, T^, Q^, T^ and r^ as functions of time. These 
equations are solved numerically by a forward difference method. The Instantaneous 
burning rate will differ from the steady-state value at a given oressure to the extent 
that Tj differs from T^. This difference can arise as a result of the changing thermal 
profile In the solid In response to the transient heating. 

3.3.3 Equations Related to the Surface Structure of Nitramlne Propellants 

Additional equations enter Into the system for nitramlne propellants because 
of the relationship of the surface structure to the burning rate (Refs. 1, 7) and 
mass flow rate. It Is assumed that the surface structure will develop In three 
stages In the course of the pressurization transient: 


a) formation of the melt layer following Ignition; 

b) disappearance of the melt layer when the break point criterion Is 
reached; 

c) formation of the surface craters following disappearance of the melt 
layer. 

The melt layer will start to form as soon at = T|yj, the melting point of 
the nitramlne. The significance of the melt layer thickness Is that It positions 
Qjjjj, the heat of fusion, which subtracts from the exothermic In the energy 
balance equations. It Is expected that the location of T|^ will move Into the solid 
with time; Tj^j Is 551 ®K for HMX, and computed steady-state surface temperatures are 
circa 1000®K. 

The second stage conmences when the critical burning rate Is exceeded, 
r^>r^j^, such that subsurface melting can no longer propagate. At this time step, 

Qy^ Is set to zero and the existing melt layer thickness Is allowed to decompose 
away In accordance with the time Integral of r^ In succeeding time steps. This 
stage Is completed when the melt layer has disappeared. 

The formation of the craters Impacts both r and In accordance with the 
steady-state model. However, the development Is allowed to proceed continuously 
with time (In the steady-state model It appears as a discontinuity with pressure). 
The changing surface area In the course of crater development Is expressed as 
follows: 




(9) 


Abp " "frozen", or cratered surface burn area 

Ajj » geometric, or Input burn area 

* normalized oxidizer surface exposed for a planar 
° M surface, from Ref. (7). 


(S /S ) c ■ normalized oxidizer surface exposed for a 
° cratered surface, from Ref. (7). 
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The parameter (Sq^/^q)p changes with time because the oxidizer (nitramlne) 
penetration depth, h, changes with time. The relationship between 
and h appears In Ref. (1). The variation of the penetration depth with time 
may be expressed as follows: 

dh « - r^) dt (10) 

'"ox * *ox <-Ws> <"> 

• A, exp (-Wj) 02) 

*^ox “ f’^tramlne monopropellant burning rate 
• binder regression rate 
Aox “ kinetics prefactor for the nitramlne 
Eox * activation energy for the nitramlne 
« kinetics prefactor for binder 
» activation energy for binder 

A limit Is placed upon the extent of penetration In accordance with tlie steady- 
state model. This limit results In a maximum (Sqjj/S^)p, or (Sqjj/Sq)|^jj. From 
the steady-state model. It Is known that the area change from 
(Sox/So)p Is the essential manifestation of the break point wherein r jwnps 
from a value Influenced by the binder to a value close to the nitramlne mono- 
propellant rate. Therefore. In the course of the break point transition. It 
may be assumed that r Is proportional to the change In This assumotlon 

Is made In order to avoid Jump discontinuities In the time-dependent model, and 
allow smooth (continuous) transition. Eq. (2). as applied during this third 
stage of surface development, becomes: 



exp 


1 




P 


r (-Ec/R) (l/T, - 

1 ^ 



^^oxV’^^ox^M 


L o a a J 


•■m 



. 


(13) 


rp> post-break steady-state burn rate 
pre-break steady-state burn rate 
(^)p ■ that value associated with Fp 
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occur at the same pressure as In steady-state. That Is, the pressures for 
r^>rjjj and r> r^^^ may be different. Thus It Is necessary to Include values of 
r which extrapolate planar surface values (ry^) to higher pressure and cratered 
surface values (Fp) to lower pressure. If the cratering actually occurs at 
some different pressure than steady-state, a proper value of r may still be used 
In accordance with the surface existing at the time. 

The additional unknowns presented by nitramlne propellants are (Sqj^)|,|, r^^^, 
rp, h, (S^^)p and A^p. Eqs. (9-12) and two equations from the steady-state 
model are Included in the forward difference scheme to provide the relationships. 

Eq. (13) substitutes for Eq. (2), and A^p Is used In Eq. (6). The changing surface 
structure Is an added means by which the Instantaneous burn rate may differ from 
the equilibrium value, and a difference between A^^ and A|^p would bear upon effective 
burning rates deduced from the closed vessel experiment. The effective burning 
rate Is obtained by multiplication of Eqs. (9) and (13). 


3.3.4 Computer Programming 

A computer program providing for the nimierical solution of the equations was 
written for operation on a UNIVAC 1108 computer. The program has been successfully 
rendered operational and applied to test cases. An Important limitation Is the 
run time which can be encountered at very high burning rates, because a rapid 
transient demands a smal 1 time step for accuracy. To minimize the Impact of this 
dependence, a scheme which automatically adjusts the time step and the solid 
propellant grid to the existing transient was successfully developed. The run 
time for a case In which the burning rate reached 20 In/sec was 15 minutes. A 
description of the computer program Is presented In Appendix B. A source deck 
and listing can be furnished upon request. 
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MODEL RESULTS 

3.4.1 Howoqaneous Propellant Test Case 

The homogeneous propellant test case selected for model validation Is "Test 
Propellant A" as specified by the JANNAF Workshop on Burning Rate Measurements, 
Memoranda of 26 January 1977 and 5 April 1977. This case was selected because of 
the extensive data provided In association with the purposes of the Workshop. 

Burning rates were reported by six different laboratories for a given set of 
specifications and closed vessel pressure-time data. However, the identity of 
the propellant was not disclosed and strand burning rates were not given. In the 
absence of strand data, the reported burning rates were used as baseline equilibrium 
data for purposes of model application. The data were reported over the pressure 
range 5Kpsi-40Kpsi . An average of the results for the six laboratories was used at 
each pressure, and an extrapolation was performed to complete the input down to 
ambient pressure. Although these data are closed vessel data, the model would 
still be useful to show the effect of the closed vessel environment, if any. If 
such effect is small, as expected, the model should reproduce these data as well 
as the measured pressure-time curve. 

The burning rate results are shown in Figure 3. The solid line represents the 
actual closed vessel data, except for pressures below 5Kpsi wherein It Is an 
extrapolation. The model prediction is represented by the circles. It is observed 
that the prediction is very close to the line between SKpsi and 40Kpsi, indicating 
a very small effect of the transient environment. In fact, the effect of the 
transient environment disappears completely for this test case above 13Kps1. 

An Interesting undulation appears in the predicted results between 4Kps1 and 
12Kps1. It is not certain whether this reflects a real mechanism, a tradeoff 
between the transient combustion zone heating and transpiration cooling during 
the pressure rise, or is simply a result of the forward difference numerical method. 
Similar undulations have appeared in computerized data reduction for closed vessel 
burn rate (e.g.. Ref. 24). 
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Figure 3. Predicted Effect of Closed Vessel Environment on Burn Rate 
of a Homogeneous Propellant 
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The most Interesting behavior appears at low pressure. The very high Initial 
burning rates are largely a consequence of the assumed Igniter. The model does 
not Include a representation of the Igniter as such, but uses an arbitrary 
Input heat flux to get the burning under way. If this heat flux remains operative 
during burning, and Is large compared to the combustion-induced heat flux, the 
burning rate will be high. The Igniter Is turned off at 1000 psi, and the 
burning rate Immediately drops because the transpiration cooling becomes excessive 
for the remaining heat Imposed. The burning rate climbs again when the transient 
heating Is re-established. Here, then. Is a second type of burning rate undulation. 

A third type (not shown here) results from pre-ignition thermal soak of the propellant, 
or dynamic "overshoot" due to a sudden Imposed compression" (25). Thus, there are 
mechanistic bases for low pressure burning rate undulations, depending upon the 
nature of the Ignition process and developing combustion-induced heating. These 
effects, although real, are undesirable for purposes of closed vessel data 
acquisition and are lumped into the area of uncertainties. As a result, burning 
rates deduced from the closed vessel generally exclude results below 5 kpsi. An 
accurate quantification of the low pressure behavior would require a model of the 
igniter. 

Predicted pressure-time results are shown together with data in Figure 4. 

Time is mt'sured from ignition time rather than zero time to remove the difference 
in ignition delay and afford a better visual comparison. Initial pressurization 
is overpredicted because of the effect of the assumed igniter. This also causes 
the predicted web to advance relative to the actual web at a given time and 
pressure. The curves then merge and the prediction becomes very good at ‘he 
higher pressures. The pressure results are consistent with the burning rate results. 

3.4.2 Nitramlne Propellant Test Case 

The nitramlne pt*ope11ant test case selected for model validation Is the 
92-m1cron HMX propellant reported by Lenschitz (Ref. 20). Much of the required 



78-6 



Figure 4. Comparison of Predicted and Experimental Pressurization for 
the Test Corresponding to Figure 3 
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Input Information Is contained within that paper. Including actual strand burning 
rate data to serve as the equilibrium baseline. The paper also reports the 
pressur^ atlon history for that test. Other Information required was furnished 
by Lenschitz. 

The burning rate results are shown In Figure 5. Plotted are the strand data 
(solid line), deduced closed vessel data (solid circles), and the predicted actual 
(open circles) and effective (x) burning rates In the closed vessel. As expected, 
there Is a significant difference between strand and closed vessel burning rates in 
the c«se of a nitramlne propellant, and the model se.^ves to explain why. 

First, let us examine the region below 2000 psi. Here, much the same be- 
havior appears as was discussed for the homogeneous propellant; the burning rate 
is first dominated by the behavior of the assumed igniter, and drops when the 
igniter Is turned off. Note, however, the appearance of a second burning rate 
undulation just below 2000 psi. This Is a consequence of a fourth mechanism for 
undulation, and will be unique to nitramlne propellants which exhibit a break 
point. What is happening here Is that a breal- ^ oint 1s being traversed, causing 
a greater burning rate than can momentarily be sustained by the combustion heating 
then In existence. The rate falls until sufficient heating develops to re-sustain 
It. If the rate fell to the extent that the surface melt layer would re-form, 
then an endotherm would reappear to drop the rate further; such a process might 
appear as a fifth mechanism for undulation. No data are reported below 2000 psi, 
but the model results are not expected to be accurate In that regime because the 
Igniter Is not modeled. 

The region between 2000 p:.1 and 10 Kpsi is Interesting because It shows how 
the predicted actual burning rates In the closed vessel differ from equilibrium 
data. Basically, the post-break-point surface structure cratering of the propellant 
Is developing continuously and gradually In the course of the transient In the 
closed vessel. This continuous development lags the step-wise cratering at the 
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Figure 5. Predicted Effect of Closed Vessel Environment on Burn Rate 
of a Nitramlne Propellant, and Comparison With Experimental 
Data 
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discrete test-to-test pressure Intervals In the strand burner. Thus, the severity 
of the break point, as measured by the exponent In this transition region, is 
lessened by the closed vessel environment. At pressures above lOKpsI, the ore- 
dieted actual burning rates In the closed vessel merge with the equilibrium strand 
data as the surface structure becomes fully developed. 

Comparison between predicted and measured burn rates in the closed vessel re- 
quires use of the predicted effective burn rate, because that parameter accounts 
for the effective propellant burn surface area due to the cratering. The predicted 
effective burn rate does not depart significantly from the predicted actual burn rate 
until pressures in excess of 3Kps1. Over the range 3Kpsi-20Kpsi , the predictions 
are generally on the low side, but agree rather well with the data. Thus, it is 
concluded that the high closed vessel burn rates at high pressure, relative to the 
strand data, are due entirely to a surface area effect and not a transient burn 
rate effect. Lenschitz had rendered a similar conclusion (Ref. 20), but thought 
that the propellant grains were breaking up. Grain breakup is not required to 
explain these results if the cratered surface structure is recognized. 

The predicted pressurization is compared with data in Figure 6. Here, the 
data were reported as rate of change of pressure versus pressure rather than as 
pressure versus time. This is a more severe test of the model because a number of 
factors enter into the rate of pressure change at a given pressure. The prediction 
falls below the data, which is consistent with the underprediction of the effective 
burn rates, but the agreement is rather good considering all of the factors involved. 

Based upon the results of this effort and the stesjy -state modeling effort, 
it would be expected that closed vessel burning rates should agree with strand 
burning rates when the '^X particle size is sufficiently fine that there is no 
break point or surface cratering. 
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Figure 6. Comparison of Predicted and Experimental Pressurization for 
the Test Corresponding to Figure 5 
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SECTION 4 

COMBUSTION INSTABILITY CHARACTERISTICS OF NITRAMINE PROPELLANTS 

4.1 BACKGROUND HYPOTHESIS 

A belief that HMX may have favorable effects upon combustion Instability 
driving, relative to AP, Is based upon the qualitative Information presented In 
Table 1, With respect to some factors which may be Important, HMX appears to 
be more similar to ammonium nitrate (AN) and potassium perchlorate (KP) than to 
AP. Propellants based upon AN or KP have not been reported to be unstable 
(Ref. 26). HMX Is similar to AN and KP In that there Is a substantial oxidizer 
melt layer and particle size effects upon burn rate are relatively weak.^ 

However, HMX surpasses AP as an energetic monopropellant. To the extent that 
the products of combustion Include particulates which may afford acoustic damping, 
HMX would tend to produce more carbonaceous residue than AN (but perhaps not 
as effective as spherical KCl droplets). At the commencement of this work, 
there was no reported systematic study of the effect of HMX on combustion driving. 

4.2 PROPELLANT SELECTION 

Standard thermochemical calculations and burning rate model calculations 
were performed to provide selection of the propellant fornwlatlons for T-burner 
combustion Instability experiments. The purpose Is to determine the effect of 
replacing ammonium perchlorate (AP) with HMX on the response function of Inert 
binder (HTPB) propellants, and of HMX addition on the response function of 
nitrocellulose (NC) propellants. In making these comparisons. It was desired 
to maintain propellant energy and burning rate reasonably constant within each 
group. Although differences can be corrected In the course of Interpreting the 
response function. It Is considered that excessive differences will render the 
comparison less meaningful due to mechanistic bias. 

T true at pressures below the HMX break point. Significant changes occur above 
the break point, which would alter the hypothesis. Pressures of Interest in 
solid rockets, as opposed to guns, are generally below these break points. 


TiMlLE 1 


PREMISES FOR COMBUSTION ST*\BILiri OF HMX PROPELLANTS 


Propensity Toward Instability Yes 

Substantial Oxldixer Melt No 

Bum Rate as Function of Strong 

Particle Size 

Condensed Phase Heat Release Strong 

Moncpropellant Plane Energy Weak 

Particulates in Combustion Zone No 




HMX 

t;o 

Ho 

7 

Yer 

Yes 

Yes 

Wepk 

Weak 

Weak 

Vet k 

None 

Strong 

Wei.k 

Non-. 

Strong 

!*■ 

Yes 

?* 


* - Depends on fuel- richness. For a given solids loading, fr. HMX propellant 
would be nore fuel-rich than an AN propellant. Although It Is possible 
to formulate such that no free carbon would be predlctet. by equlllbrltnn 
thencochealstry, residues are sometimes observed. 
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AP and HMX produce different thermochemical effects In HTPB binder because 
AP Is an oxidizer whereas HMX Is an energetic monopropellant. HMX propellants 
tend to be low burning rate propellants* but significant burning rate control Is 
available In AP propellants through particle size adjustment. Tradeoff study 
produced the two groups of selected HTPB propellants listed In Table 2: 
a low burning rate group (XA) and a high burning rate group (XB). It Is difficult 
to match AP propellant burning rates with HMX alone, so the burning rate of 
formulation XA-3 Is lower than desired for Its group although Its energy level 
Is adequate. The two mixed oxidizer propellants are more closely related to their 
AP analogs In burning rate and energy. 

It was found that the addition of meaningful quantities of HMX to an inert- 
plasticized NC binder produced considerable Increases In propellant energy. Use 
of a relatively stable energetic plasticizer (TMETN) reduced the impact of 
HMX addition on propellant energy, and was considered more acceptable for 
research purposes than nitroglycerine. The resulting propellants comprise the 
third group (XC) listed in Table 2. These propellants were not tested in the 
course of this program because additional testing of the HTPB propellants was 
found to be necessary In order to acquire sufficient meaningful data. 

4.3 EXPERIMENTAL RESULTS 

Stability tests were carried out In 2 l/2-1nch diameter (I.D.) T-burners 
coupled directly to a surge tank. The propellant configuration used was a 3/8- 
Inch thick disk. Data were obtained using the growth-decay technique (Ref. 27). 

For each propellant, tests were conducted at nominal frequencies of 500, 900, 
and 1900 Hz and at pressures of 500 and 1000 Psi. Several tests were performed 
at each condition. Many of the tests did not produce oscillations, or the data 
did not conform with T-burner standards, and those tests are excluded from the results. 
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TABLE 2 


PROPELLANT FORMULATIONS FOR T-BURNER TESTS 



XA-1 

XA-2 

XA-3 

XB-^ 

XB-5 

XC-6 

XC-7 

AP(9ji) 

- 

- 

- 

41.25 

31.25 

- 

- 

AP(50y) 

40.00 

20.00 

- 

- 

- 

- 

- 

AP(90y) 

- 

- 

- 

41.25 

31.25 

- 

- 

AP(200y) 

40.00 

20.00 

- 

- 

- 

- 

- 

HMX(4y) 

- 

20.00 

40.00 

- 

22.50 

- 

50.00 

HMX(15y) 

- 

20.00 

40.00 

- 

- 

- 

- 

R-45/HTPB 

20.00 

20.0c 

20.00 

17.50 

15.00 

- 

- 

NC/TMETN 

- 

- 

- 

- 

- 

100.00 

50. CO 

(T„/M)^^^ 

£ 

10.34 

10.30 

10.32 

10.70 

10.77 

10.49 

10.86 









■'1000'^"/''“’ 

0.26^ 

0.21^ 

0.16^ 

0.48^ 

0.43^ 

0.30 

0.33 


1 - Measured Values 
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The results for the low rate and high rate HTPB propellants are given In 
Figures 7 and 8, respectively. The results are plotted as the sum of the measured 
growth and decay coefficients vs. frequency. The propellants were all relatively 
stable In that the measured growth coefficients exceeded the decay coefficients 
only for the two high rate propellants at the 1900 Hz, 500 ps1 test condition. 

Such behavior may subject the absolute results to criticism, but the relative 
results are considered valid. 

The low rate XA-1 propellant was the weaker acoustic driver of the two 
propellants which did not contain HMX. The low rate results from the use of 
coarser particle sizes. Its results appear to be relatively Independent of 
pressure. The minimum between 500 Hz and 1800 Hz Is considered real because only 
one test out of seven conducted at 900 Hz produced measurable acoustic driving, 
whereas most of the excluded data at the other frequencies stemmed from poor 
raw data. 

The two low rate propellants containing W1X exhibited unusual results. At 
500 Hz, the pressure traces showed an Immediate growth of oscillations following 
Ignition to amplitudes of about 2% of the mean and then a slow decay over the 
duration of the burn. At 1900 Hz, the tests exhibited no oscillations at all. 

At 900 Hz, the 500 psi tests were like those at 500 Hz and the 1000 psi tests 
were like those at 1900 Hz. 

Post-test examination of the XA-3 tests showed the T-burner to be full of a 
layered, carbonaceous char. In the 1900 Hz tests, this char appeardd to have 
been pulverized. For the lower HMX concentration XA-2 propellant, the char was 
less extensive, merely coating the wall of the T-burners rather than filling the 
Interior. The cha» s presumed to be due to the fuel -rich nature of the two low 
rate propellants containing HMX. It was somewhat surprising because prior 
experience with similar propellants did not encounter char. A review of the 
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FREQUENCY, Hz 


Figure 8. T-Burner Results for High Rate 
HTPB Propellants 
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formulations revealed that the prior propellants were much more heavily plasticized 
gun propellants, whereas the propellants formulated here were more consistent 
with the structural integrity requirements of rockets. Another consequence of 
the reduced plasticizer, or the char, was that the burning rate was about 20% 
lower than expected. The char Is Indicative of non-equilibrium thermochemistry 
In that there should be no free carbon at a solids loading of 80% HMX. In any 
event, the contribution of HMX to stability could not be identified between 
combustion driving or mechanical suppression from these tests. 

The high rate propellant data in Fig. 8 is irore conclusive as to the role 
of HMX. These propellants did not produce char. Both XB-4 and XB-5 showed 
clear pressure and frequency effects, increasing in d»*iving strength with de- 
creasing pressure and increasing frequency. The XB-5 HMX-containing propellant 
would not drive oscillations at all at the 1000 psi - 500 Hz test condition, 
and its data at 500 psi and 1000 psi fall below the respective data for the 
all-AP oxidizer XB-4 propellant. The measured damping coefficients at each 
frequency were about equal for the two propellants, indicating that the differences 
in acoustic driving are real. More data are required for clean-burning oropel- 
lants, and testing of the XC group is included in future plans. 
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SECTION 5 

EFFECT OF NITRAHINE INGREDIENTS ON STEADY-STATE BURNING RATE 
(Data Supplement to Annual Progress Report AFOSR-TR-76-1136) 

5.1 PREFACE 

Strand burning rate data supplementing the data reported In the FY 1976 
Annual Progress Report (Ref. 1) were furnished by AFATL following the publica- 
tion of that Report. These data are presented and discussed In the following 
subsections. It Is concluded that the data provide further confirmation of the 
results and predictions of the steady-state model developed In the course of 
that research. 

5.2 EFFECT OF HMX PARTICLE SIZE IN ACTIVE BINDER PROPELUNTS 

Figure 5 of Ref. (1) showed the effect of HMX particle size on the burning 
rate of nitroplasticized polyurethane (NPPU) propellant. Analogous data for NC 
propellant are now available and are presented In Fig. 9. The same basic trends 
appear as In the NPPU propellant, as was expected. The fine propellant exhibits 
a relatively uniform pressure exponent over two orders of magnitude of pressure. 

The coarse propellant begins to depart significantly at a pressure above IKpsI, 
with an upward break In pressure exponent covering the range 1.3Kps1 to 2.5Kos1. 

A second upward break appears to begin at 4Kps1 and continue to high pressure. 

As a result, the burning rates of the coarse propellant have a muUiole zig-zag 
appearance and are higher at high pressure. The burning rates of the NC propellants 
are higher than those of their NPPU analogs, which Is attributed to the higher 
energy and burning rate of the NC binder. 

A similar trend Is shown by the mixed HMX-TAGN propellants In Fig. 10. A 
portion of the data for the propellant containing the fine HMX appeared In Fig. 7 
of Ref. (1). No data were available for the propellant containing the coarse 
HMX. Results for both propellants are now complete. The propellant 
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Figure 10 Effect of HMX Particle Lise In Mixed HMX>TAGN, NC Propellants 
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containing fine (MX exhibits a relatively uniform shape over a wide pressure 
range. The propellant containing coarse HMX exhibits an upward break In pres- 
sure exponent covering the range IKpsI to 3Kps1, such that Its burning rates are 
higher at high pressures, as was expected. Comparing Figs. 9 and 10, the effect 
of the TAGN Is to raise burning rates, more so at low pressure than at high 
pressure, such that the overall net pressure exponent Is reduced. This effect 
has been attributed to the faster kinetics of TAGN decomposition and combustion 
as a monopropellant (Ref. 1). 

Figs. 9 and 10 Illustrate the effectiveness of fine particle size as an 
approach to achieve desirable slope uniformity. Fig. 9 confirms that it Is 
possible to do so with HMX alone If the particle size Is sufficiently fine. 

Fig. 10 confirms that TAGN addition Is helpful provided that the HMX component 
does not exhibit a slope break. 

5.3 EFFECT OF TAGN IN ACTIVE BINDER PROPELLANTS 

Data presented in Figs. 9 and 10 allow the completion of Fig. 7 of Ref (1). 
This completed figure Is presented as Fig. 11. The effect of TAGN substitution 
for HMX, In Increasing amounts. Is thereby Illustrated. It was predicted that 
the data for the 29% TAGN propellant would Interweave with the data for the 
50% TAGN propellant at the higher pressures. This prediction Is verified by 
the data. A break point appears for the 50% TAGN propellant at lOKpsl. This 
was explained by the burning rates reaching such a level as to trigger the break 
In the HMX component, and was described as a disadvantage of excessive TAGN 
addition. In view of that explanation, the 29% TAGN propellant also should 
exhibit a break point when the burning rates reach the same level. This Is not 
so obvious In Its high pressure data, perhaps because the region Is not as well 
characterized with data points. In any event. It Is confirmed that. TAGN addition 
raises burning rates, and that optimum amounts are helpful as an approach to 
achieve desired slope uniformity. 
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Figure 11 Effect of TAGN in NC Propellants 
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There has been a question as to whether or not propellants containing TAGN 
alone would exhibit a break point. The uppermost data In Fig. 11 are for such 
a TAGN propellant. The model predicted that a break point exists near the region 
where considerable data scatter appears, and the scatter was Interpreted to be a 
manifestation of the break. The model also predicted that the second upward 
break In the 50% TAGN propellant is due to the TAGN. The question, however, 
could not be resolved from the data. Complete data are now available for the 
63% TAGN, NPPU propellant (cf. Fig. 6 of Ref [1]). These data are shown to- 
gether with the Fig. 11 data in Fig. 12. A break point for this TAGN propellant 
clearly appears at lOKpsi. Therefore, It can be concluded that TAGN propellants 
can exhibit break points. Ironically, this break point can be mitigated by 
lowering the burning rates, through H|^ substitution, which would yield the GAU-8 
propellant. The shape of the NPPU propellant burning rate curve lends credence 
to the interpretation of a high slope (break) region In the case of the NC 
propellant. The NC propellant has higher burning rates because of the greater 
TAGN concentration and the more energetic binder. 

5.4 MATCHING NITRAMINE AND ACTIVE BINDER 

It was pointed out In Ref. (1) that matching the burning rates of the active 
binder and nitramlne as monopropellants could yield a propellant having a uniform 
pressure exponent. The rationale was that a break point would have no way to 
manifest Itself, and so would not appear In the burning rate data. In support 
of this hypothesis, data and results of model calculations were presented In 
Figs. 31-33 of that Report. These figures showed trends as the nitramlne 
burning rate approaches and falls below the binder burning rate.. It was suggested 
that a propellant combining fine EDNA and a tailored active binder having the 
same burning rates as EDNA would be a practical test of the approach. 
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Figure 12 Burning Rates of TAGN Active Binder Propellants 
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Unfortunately* the propellant that was furnished combined coarse EDNA and 
the baseline NC binder. As a result* the binder burning rates are lower than 
the EDNA burning ratis In accordance with the Figure 31* and the break point 
mechanism would be triggered at a low pressure. With this propellant* however* 

It Is possible to complete the EDNA propellant portion of Fig. 31 and these data 
are shown In Fig. 13. The premise that the post-break propellant burning rates 
follow the nitramlne monopropellant burning rates Is confirmed. Without data at 
lower pressures, where a break would be expected* It does not appear to matter 
that the binder burning rates are low. A transition or jump upward to the EDNA 
line Is not seen over this pressure range* such that this EDNA/NC propellant 
appears to achieve a desirable result. 

An illustration of binder tailorablllty to Implement the foregoing approach 
is presented in Fig. 14. The lower dashed line is for nitrocellulose (NC) 
containing an inert plasticizer. The upper dashed line is for nitrocellulose 
containing a very energetic plasticizer. The line fitting AFATL data is for 
nitrocellulose containing an energetic (intermediate) plasticizer. The range of 
burning rates is observed to cover a factor of 3. A larger range would be available 
by adjusting ingredient proportions* subject to energy and processing limitations; 
a practical range of 5 is indicated by existing active binders. The intermediate 
binder shown in Fig. 14 has burning rates quite close to EDNA burning rates, and 
would be an attractive candidate. Note, however* that the details of the binder 
data indicate some break point mechanism in the binder which might be reflected 
by the propellant using it. It is difficult to coinnent upon break points in active 
binders in general because the data are usually reported as straight lines. 
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Figure 13. Comparison of Monopropellant Burning Rates and EDNA 
Propellant Burning Rates 
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SECTION 6 

PRESENTATIONS AND PUBLICATIONS 


The following presentations and publications have been generated by this 
research contract: 

(1) Cohen, N.S. and Strand, L.D., "Nitramine Propellant 
Research", 13th JANNAF Combustion Meeting (CPIA Pub- 
lication 281, Vol. I, Dec. 1976) pp. 75-87. 

(2) Cohen, N.S. and Strand, L.D., "Nitramine Smokeless 
Propellant Research", 1977 AFOSR/AFRPL Rocket 
Propulsion Research Meeting, Lancaster, CA (Mar, 1977). 

(3) Cohen, N.S., Strand, L.D., and Price, C.F., "Analytical 
Model of the Combustion of Multicomponent Solid Propellants", 
AIAA Paper 77-927, AIAA/SAE 13th Propulsion Conference, 
Orlando, FL (July, 1977). 

(4) Cohen, N.S. and Strand, L.D., "Transient Combustion 
Response Characteristics of Nitramine Smokeless 
Propellants", 14th JANNAF Combustion Meeting , Colorado 
Springs, CO (Aug., 1977), publication pending. 
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SECTION 7 
CONCLUSIONS 

A comparison of strand and closed vessel burning rate data for homogeneous 
propellants reveals that effects of pressurization are minor, and that the two 
methods should give consistent results at pressures following the Ignition 
transient In the closed vessel. A description of burning rates during the 
Ignition transient In a closed vessel would require a model of the Igniter. 

The burning rates of nitramlne propellants as measured In a dosed vessel 
win differ from strand data for those propellants which exhibit a break point. 

The effect of the transient environment In a closed vessel Is to create a time 
lag In the development of the post-break-point cratered surface structure, which 
results In a less severe exponent shift than shown by the strand data. The 
reduced severity of the exponent shift will have a practical benefit In transient 
devices. However, the cratered surface structure will maintain high mass burning 
rates and exponents because of the increased burn surface area. Methods to avoid 
or mitigate the cratering, as discussed In the Annual Progress Report for FY 1976 
(Ref. 1), are still recommended. Additional data have been presented here which 
confirm the conclusions of that Report. 

Several mechanisms have been Identified which might explain burning rate 
undulations often observed from closed vessel data reduction. In particular at 
low pressure: Ignition overshoot, termination of the Igniter, tradeoff between 

developing combustion zone-induced heating and Increasing surface transpiration 
plus. In the case of nitramlne propellants, break-point overshoot and melt layer 
reformation. 

HMX-contalnIng composite propellants are more stable than comparable AP 
propellants. One mechanism may be related to mechanical suppression due to the 
generation of solid material combustion products In certain fuel-rich formulations. 
However, It also Is true that HMX reduces the combustion driving In clean-burning 
formulations tested thus far. Additional study Is reconmended. 
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APPENDIX A 
NOMENCLATURE 








t 


kinetics prefactor for propellant decomposition 
propellant burn surface area 

nitramlne propellant burn surface area following cratered surface development 

kinetics prefactor for binder decomposition 

kinetics prefactor for nitramlne decomposition 

surface area of closed vessel wall exposed to heating 

propellant heat capacity 

gas heat capacity 

activation energy for propellant decomposition 
activation energy for binder decomposition 
activation energy for nitramlne decomposition 
Igniter heat flux 

depth of surface cratering In nitramlne propellants 

heat transfer coefficient for wall heatlno 

propellant thermal conductivity 

mass of gas In chamber 

burning rate pressure exponent 

pressure 

rate of change of pressure 

heat flux at propellant surface 

propellant heat of decomposition 

rate of heat loss to closed vessel wall 

steady-state burning rate at P 

regression rate of binder 

steady-state post break point burn rate at P 

Instantaneous propellant burn rate 

steady-state pre-break point burn rate at P 

regression rate of nitramlne 

gas constant 

universal gas constant 

fraction of exposed nitramlne surface on a planar melt surface 
fraction of exposed nitramlne surface on a cratered surface 

value of associated with Fe 
ox r 

time 


r- .‘-it,- -iMt 
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Tp propellant adiabatic flame temperature 
Tg Instantaneous gas temperature 

Tp Initial (ambient) temperature 

Instantaneous propellant surface temperature 
T steady-state propellant surface temperature associated with r 

closed vessel wall temperature 
V Instantaneous gas volume In the closed vessel 

w propellant web 

X distance from propellant surface Into propellant 

e radiation emissivlty of gas 

n covolume 

y ratio of specific heats 

K propellant thermal diffusivlty 

p propellant density 

(j- Stefan-Boltzmann constant 

i steady-state flame height associated with r at P (dimensionless) 
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APPENDIX B 

COMPUTER PROGRAM DESCRIPTION 


B-1 SUBROUTINES 

The computer program file is identified as KT, and its map of linked 
• elements for execution is identified as KRTT. The component subroutines 

are listed as follows: 


Segment 

MAIN PROGRAM 

Subroutine 

EVES 

Subroutine 

SETUP 

Subroutine 

DIFEQ 

Subroutine 

TSWALL 

Subroutine 

TSPROP 

Subroutine 

PRINT 

Subroutine 

FIND 


A schematic of the computer program elements is shown in Fig. B-1. 

Segment MAIN is the master control element. This element reads input, 
computes certain parameters of an input nature, writes the input, and 
then transfer control to EVES. Subroutine EVES is a simultaneous 
differential equation solver that uses the Adams preoictor-corrector 
method. EVES sets the time step, predicts ahead, compares the results 
emd certifies accepted answers. Subroutine SETUP initializes all of the 
variables, and is the first subroutine called by EVES. Thus SETUP provides 
the initial conditions for the time-dependent problem. Subroutine PIFEO 
computes the time derivatives of qas pressure, gas volume, gas weight, 
gas temperature and distance burned. DIFEQ contains the time -dependent 
conservation equations which compute these differentials. EVES calls 
DIFEQ after the problem has been initialized. DIFEQ calls TSPROP to 
obtain the propellant combustion contributions to the conservation 

, 
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a) MAX. RUN TIME EXCEEDED 
OR b) MAX. PRESSURE EXCEEDED 
OR c) WEB BURN-OUT ACHIEVED 


Figure B-1. Schematic of Computer Program Elements 
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equations, and TSWALL to obtain the heat loss to the wall which enters 
into those equations. TSPROP computes the time-dependent temperature 
profile in the solid propellemt, the burning rate which is a function 
of surface temperature, and includes mechanisms by which nitramine behavior can 
influence the burning rate. A finite difference approximation of the 
transient heat equation is used to compute the temperature profile, which 
develops in response to igniter-induced heating and then to combustion 
heating. TSPROP may control the time step when high burn rate leads to 
steep temperature profiles. TSWALL ccxnputes the convective and radiation heating 
j of the wall, and the temperature profile in the wall, in determining the 

heat loss. TSPROP and TSWALL return answers to DIFEO, and DIFEQ returns 
answers to EVES. EVES then calls PRINT to update the accepted answers and to 
write them at pre-selected time intervals. Subroutine FIND is an auxilicury 
subroutine that is used to interpolate. 

The computations involving TSPROP, which is the essence of the 
combustion model, warrant a more detailed discussion which is provided in 
Section 2 of this Appendix. A source deck and listing of the ccxnputer 
program can be furnished upon request; it is liberally embellished with 
comment cards to assist the user. A manual for inputs is provided as 
Section 3 of this Appendix. Output is discussed in Section 4. 

B-2 THE COMBUSTION MODEL CALCULATIONS (TSPROP) 

The following sequential discussion describes the computations performed 
in TSPROP. 

1. The time step is first Jeteruined by EVES. 

I 

2. Tabul£u: strand burning rates are initialized for the beginning of 
this time step. 
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3. The propellant web calculational mesh and counter are updated 
for the web which was burned away in the previous time step. 

4. A further euljustment to the mesh is made depending upon whether 
the reference mesh is greater than 1/10 the thermal wave thickness 
of the previous time step, or less than 1/25 of that thickness. 

Too coarse a mesh impairs accuracy, too fine a mesh will increase 
computer time by requiring a smaller time step. If the mesh is 
adjusted, the temperature profile is redefined by logarithmic 
interpolation to conform to the new mesh. This technique of 
adjusting the mesh has provided considerable savings in computer time. 

5. The time step required by the mesh is compared to the time step 
required by EVES, and the smaller value controls. 

6. The igniter -induced heat flux is gradually diminished to zero when 
sufficient burning has been achie’ ed. The method of representing 
the igniter is arbitrary. There is no model of the igniter as 
such, only an input heat flux to start. 

7. The heat flux to the propellant is computed, consisting of igniter 
and combustion contributions. The actual heat flux is referenced 
from that which would exist under equilibrium burning at the pres- 
sure of the present time step, and therefore assumes a quasi-steady 
gas combustion zone. The appropriate tabular strand burning rate is 
updated to the present time step for this purpose. If the propellant 
is a nitramine propellant, and the actual burning rate of the 
previous time step is below the breakpoint burning rate, then the 
exothermic surface heat release term is modified by the heat of 


fusion of the nitramine. 
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8. Using a finite-difference expression, the resulting heat flux 
is converted into an 2 urtificial surface temperature that is 
used to drive the thermal profile in the depth of the solid. 

9. The finite difference form of the transient heat equation is used 
to compute the temperature at each mesh location in the depth of 
the solid, including the true surface temperature. The rate of 
change of temperature is computed and multiplied by the time step. 

10. The actual burning rate is computed from the surface temperature 
using an Arrhenius expression. 

11. Prop>ellant density is known from inputs. The present burn eurea 
is known from subroutine DIFEQ. Knowing the actual burn rate 
the propellant weight flow rate is then computed. 

12. If the propellant is hcxnogeneous , this burn rate and weight 

flow rate are returned to DIFEQ at this point. If the propellant 
is a nitramine propellant, additional computations are performed. 

13. A test is made comp^u:ing the actual burning rate with the break -point 
burning rate. If the actual burning rate is less, the melt layer 
thickness of the nitramine is determined by the location of the 

melt temperature in the thermal profile. If the actual burning 
rate Is now less but was at one time greater than the break point 
rate, such that the surface is now cratered, then an additional 
computation is performed to begin or continue to knit the surface 
back to a planar melt. This latter condition can result from de- 
pressurization or reductions in heat flux. 

14. After computing the melt layer thickness, the program returns to 
DIFEQ if the surface la now or is once again, a planar melt. If the 
surface has not yet knitted completely, then adjustments are made 
according to the extent of cratering as discussed below. 
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15. If the actual burning rate is greater than the break-point 
burning rate in step 13, then the melt layer thickness is 
reduced by tha product of the burning rate and the time step. 

If the melt layer thickness is still greater than zero, the progreun 
returns to DIFEQ. If the melt layer has disappeared by this 
mech 2 mism, then the extent of cratering is conoputed according 
to the difference between the nitreunine and binder burning rates. 
These ingredient burning rates are calculated from Arrhenius 
expressions for the present surface temperature. If the binder 
is an active binder, its rate is determined from the updated 
tabular input. 

'* 16. If the surface is cratered, either as a result of step 15 or 

step 14, the actual burning rate is upgraded in accordance with 
the mechanism contained in the steady-state model. The 
parcuneters which produce this upgrade are computed from the 
known extent of cratering and tabular burn rates as presently 
updated. Associated with this upgraded burn rate is an upgraded 
burn area, also computed from the cratering. Thus em upgraded 
weight flow rate is computed, and these upgraded values are 
returned to DIFEQ instead of the values from Step 12. 

B-3 CARD INPUT SEQUENCE AND STANDARD VALUES 

The following describes input cards in sequence, with appropriate 
comments and numerical values which have been used. 
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JOB CARD 

SEQUENCE OF CONTROL CARDS 
EXECUTION CARD 


CARD 

COLUMN 

ENTRY 

SYMBOL 

UNITS 

1 

2 

Type of Propellant 

ICASE 

- 



Canment; ICASE = 1 for homogeneous 

propellants 




ICASE •= 2 for nitramine propellants 


2 

3 

Title of Case or Group 

ITIT 


3 

2 

Initial Chamber Gas Volume 

VINIT 

. J 

in 2 


12 

Initial Pressure 

PINIT 

Ib/in 


22 

Initial Temperature 

TINIT 

°K 2 


32 

Igniter Heat Flux 

FEXT 

cal/cm -sec 



Comment: Values of 100-200 have been used. 



42 

Reference Calculational Grid 





size of Propellant 

DELXO 

in. 



Comment: A value of 0.0002 





has been used. 




52 

Calculational Grid Size for Wall 

DELR 

in. 



Comment: Divide wall thickness 





by 10 



3 

62 

Chamber Diameter 

DIA 

in. 

4 

2 

Surface Area of Chamber Wall 

AWALL 

"5* 


12 

Thermal Diffusivity of Wall 

ALPHAN 

cm /sec. 


22 

Flaune Temperature of Propellant 





at 1000 psi 

TFP 

°K 



Comment: Divide isochoric flame 





temperature by GAMMA 




32 

Gas Molecular Weight 

MW 

gm/mol 


42 

Ratio of specific heats 

GAMMA 

- 


52 

Specific heat of gas at 10,000 
psi 

CGP 

cal 

o— 

gm- K 


62 

Specific heat of solid propellant 

CSP 

cal 

5 

2 

Covolume 

COVOL 

gm-°K 
in /lb. 


12 

Weight Fraction of Nitramine 

ALFA 

- 


22 

Mean Particle Size of Nitramine 

DAVE 

Microns 


32 

Heocimum Surface Penetration 
Associated with Post-Break -Point 
Surface Cratering 

HDMAX 




Comment: Use a value of 3. for sizes greater than 

100 microns 


Use a value of 1. for sizes less than 10 microns 
Use a value of 2. for intermediate sizes ^ 

42 Nitramine density RHOX gm/cm^ 

52 Binder density RHOF ^qm/cm 

62 Propellant Thermal Diffusivity ALPHAP in /sec. 
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CARD 

6 


COLUMN 


ENTRY 


SYMBOL UNITS 


2 


12 


22 


32 


42 


52 


62 


2 


12 


22 

32 

42 


Reference Propellant Surface Temperature 

Temperature TSREF °K 

Cooment: For hcxnogeneous propellants, select 

a value consistent with the burning rate 
at 1000 psi and the kinetics expression 
for surface decomposition} an example is 
523 K for 0.35 in/sec. and an activation 
energy of 50000 cal/mol. Surface temperatures 
for homogeneous propellants are in this 
vicinity. For nitramine propellcuits, do 
the same but with a pre-break-point burning 
rate; for example, 926°K for 0.13 in/sec and 
50000 cal/mol. Nitrcutiine surface temperatures 
au:e circa 900°K, 2 md may be acquired from the 
steady-state model. 

Refe^sncc burning rate associated with rref in/sec. 

TSREF 

Comment: See comment for TSREF 

Nitramine Melting Point TMELT °K 

Comment: The value for HMX is 551. 

The value for RDX is 467. 


Activation Energy for Nitreunine or Hcwno- 
geneous Propellant Decomposition ESP 

Comment: A value of 50000. has been 

used for eacn. 

Activation Energy for Binder 
Decomposition EF 

Comment: Use a value of 16900. for 

HTPB. For active binders, 
input O. 

Prefactor for Nitramine Decomposition ASP 

Comment: A value of 0.5000 E -»■ 10 has 
been used for HMX and RDX. 

Prefactor for binder decomposition AF 

Comment: A value of 299 has been used for 

HTPB. For active binders, input 0. 

Surface Heat Release QSURP 

Comment; A value of -225. has been 
used to represent HMX emd energetic 
homogeneous propellants. 

Nitramine Heat of Fusion OLM 

Comment: A value of 132.5 is used 

for HMX and RDX. 

Nitramine Thermal Diffuslvitv XKAPO 

Comments: A value of 0.0011 has 

been used. 

Binder Thermal Diffusivity XKAPF 

Comment: A value of 0.0011 has 

been used. 

Nitramine Thermal Conductivity XLAMO 

Comment: A value of 0.00049 has been 

used. 


cal/mol 

cal/mol 

_3E*2_ 
cm -sec 

gm 

2 

cm -sec. 
cal/gm. 

cal/gm 

2 , 

cm /sec. 

cm / sec . 
cal 

o 

cm-sec- K 
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CARD 
7 Cont'd 


9 


NEXT 

NEXT 


NEXT 

NEXT 


NEXT 


COLUMN 


ENTRY 


SYMBOL 


52 Binder Thermal Conductivity 

Comment; A value of 0.00044 
has been used. 

62 Wall Thermal Conductivity 

2 Number of Points in the Table of 

Burn Area vs Web 

2 Web Burned 

Coiranent: The first table entry 
should be zero 
12 Burn eurea 

Comment: This is the total bvurn 
area at zero web burned. 

REPEAT FOR NWEB CARDS TO DEFINE AREA VS WEB. 

2 Number of Points in the Table of 

Pre-Break-Point Burning Rate 

vs. Pressure 

2 Pressure 


XLAMF 


XKN 

NWEB 

WEBIX 

ABX 


NPl 


PMX 


12 Pre-Break Point Strand Burn Rate RMX 


Comment: Cover the full pressure 
range with the series of caurds, 
extrapolating to PINIT and to the 
maximum pressure expected. The 
series may be omitted for homogeneous 
propellants, whence the NPl card 
should be omitted also. 


REPEAT FOR NPl CARDS TO DEFINE PRE-BREAK-POINT 
BURNING RATE VS PRESSURE 


2 Number of Points in the Table of 

Post-Break-Point Burning Rate vs 
Pressure NP2 

2 Pressure PFX 

12 Post-BreaUc-Point Strand Burn Rate RFX 

Comment: Cover the full pressure 
range with this series of cards, 
extrapolating to PINIT and to the 
maximum pressure expected. This 
series may be omitted for homogeneous 
propellants, whence the NP2 card 
should be omitted also. 


REPEAT FOR NP2 CARDS TO DEFINE POST-BREAK-POINT 
BURNING RATE VS PRESSURE 

2 Number of Points in the Table of 

Active Binder Burning Rate vs 
Pressure NP3 


UNITS 

cal 

o 

cm-sec- K 

cal 

cm-sec- K 


in 


in 


2 


Ib/in. ^ 
in/sec . 


Ib/in.^ 

in/sec 


REPBODUCroiLITY OF THE 
ORIGINAL PAGE IS POOR 
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CARD 


COLUMN 

2 

12 


ENTRY SYMBOL UNITS 

Pressure PBX Ib/in^ 

Active Binder Strand Burn Rate RBX in/sec. 

Comment: Cover the full pressure range with this 
series of cards, extrapolating to PINIT and to 
the oiaximum pressure expected. This series may 
be omitted for inert-binder propellants, but is 
a must for homogeneous propellants. 


REPEAT FOR NP3 CARDS TO DEFINE ACTIVE BINDER BURNING RATE VS PRESSURE 


LAST CARD 
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B-4 OUTPUT 

The initial portion of the output displays the input as groups labeled 
computation parameters, gas parameters, propellant parameters and the tabular 
input peurameters. The answers are displayed in columns consisting of time, 
web burned, pressure, burn rate, propellcuit surface temperature, melt layer 
thickness, crater penetration depth, fraction of exposed nitramine surface area, 
burn area and effective burn rate. Except for burn area, the last five parameters 
are associated with nitramine propellants only. The effective burn rate is the 
actual burn rate multiplied by the ratio of the cratered surface burn area to 
the ideal geometric burn area, and would be the value of burn rate deduced from 
a closed vessel test in the absence of knowledge of the actual (cratered) burn 
eurea. The burn area which is output is the geometric (uncratered) burn area. 

Answers are displayed in time increments of 0.1 msec. 

Computer run time may be controlled by a time stop or max time in 
accordance with particular facility usage, and is usually included in a control card. 
The run will automatically stop when the web burns out or the pressure of 40000 psi 
is exceeded. This maximum pressure appears in subroutine PRINT, and could be 
changed if desired. Full runs have ranged from 5 min. -15 min. CPU on a 
UNIVAC 1108. 


m 


